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Abstract. We study superconductor-ferromagnet bi-layers, not only for s-wave 
but also for d-wave superconductors. We observe oscillations of the critical 
temperature when varying the thickness of the ferromagnetic layer for both 
s-wave and d-wave superconductors. However, for a rotated d-wave order 
parameter the critical temperature differs considerably from that for the unrotated 
case. In addition we calculate the density of states for different thicknesses 
of the ferromagnetic layer; the results reflect the oscillatory behaviour of the 
superconducting correlations. 



PACS numbers: 74.62.Yb,74.78.Fk,74.81.-g 
1. Introduction 

The interplay of superconductivity and ferromagnetism has been studied for many 
years E]. Nowadays the focus is on hybrid-structures of superconductors and 
various ferromagnetic materials. For bi-layers consisting of a superconductor and 
a ferromagnetic metal an oscillation of the critical temperature has been found - 
experimentally |S] as well as theoretically 00- when increasing the thickness 
of the ferromagnetic layer. Similar observations have been made for superconductor- 
ferromagnet multi-layers [HI Ej- For two superconductors which are coupled via a 
thin layer of a ferromagnetic metal an oscillation of the critical current when varying 
the layer thickness has also been reported [HIE]; depending on the layer thickness 
such systems can be 7r-junctions. Recently a layer of a ferromagnetic metal attached 
to a bulk superconductor has been considered, resulting in the density of states 
showing oscillations ^2 El El H4| . too. The origin of these effects is a state 
in the ferromagnetic metal which is similar to that proposed for a ferromagnetic 
superconductor by Larkin and Ovchinnikov and Fulde and Ferrell pQ: in the 
presence of a spin exchange field the superconducting order parameter is spatially 
oscillating. A LOFF-like state can be present in a ferromagnet where superconducting 
correlations enter via the proximity effect 0. 

Up to now only s-wave superconductors have been studied. In view of possible 
applications, it is also important to consider an order parameter with <i-wave 
symmetry, which presumably is realized in the cuprates; there the crystallographic 
orientation fixes the direction of the <i-wave lobes. For a related experimental 
investigation see |T5] . 
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In this work we extend the theory to d-wave superconductors, where we have to 
account for the anisotropy of the order parameter. This means that in the quasi- 
classical framework, which we use, the Usadel equation is no longer applicable, and 
we have to use the Eilenberger equation instead. It should be mentioned that the 
Eilenberger equation is also needed for the description of s-wave superconductors in 
the clean limit, which was pointed out for the current problem in |16| . 

In the following we study a bi-layer of a superconductor and a ferromagnetic metal 
as presented in figure ^ As the structures we are interested in are three-dimensional 
our mean field approach is of sufficient accuracy. We consider the behaviour of the 
critical temperature and of the density of states. In particular we compare the results 
for different order parameter symmetries, namely of the s-wave and ci-wave type. First 
we briefly introduce the fundamental quasi-classical equations. Afterwards we present 
the results, and we finish with a short conclusion. 

2. Method 

To study superconductors in the vicinity of boundaries we apply the theory of quasi- 
classical Green's functions in thermal equilibrium |17l I18j . The Green's functions are 
determined from the Eilenberger equation 



Here f, is the direct product of the i Pauli-matrix in Nambu space and the identity 
in spin space; vice versa &i is the direct product of the i th Pauli-matrix in spin space 
and the identity in Nambu space. Consequently the Green's function, g, has a 4 x 4 
matrix structure. For our purpose it is sufficient to choose the orientation of the 
internal spin-exchange field of the ferromagnetic metal in z-direction, J(r) = I(r)e z , 
which leads to the term I(r)a 3 in the Eilenberger equation. 

The superconducting order parameter, which we assume to be spin-singlet, reads 



[t 3 E + I(r)a 3 + iA(p F , r) - E(r), g(E, p F ; r)] + 



+iv F ■ d r g(E, p F ;r) = . 
and must fulfil the normalisation condition 



(1) 



[g(E,p F ;r)} 2 = 1 . 



(2) 




(3) 



with 



A(p F ,r) 



r?(p F )A(r). 



(4) 




ds 







Figure 1. A bi-layer consisting of a superconductor of thickness dg and a 
ferromagnetic metal of thickness dp. The interface between both materials (at 
x = 0) is assumed to be completely transparent, and the sample is translationally 
invariant in y- and z-direction. 
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The symmetry of the order parameter is determined by the basis function, r](pp): 

{1 s-wave(f? s ) 
(Pfx ~ Pf v )/Pf unrotated d- wave (rj d ) . (5) 
PFxPFy/PF 45°-rotated d-wave (rj' d ) 

The order parameter must be determined self-consistently via 

A(r) = -7rW T E f Tr - 

\E n \<E c 

where (. . .) p , denotes an average over the Fermi surface, which is assumed to be 
spherical in the s-wave case and cylindrical in the d-wave case (this is justified for the 
layered cuprate superconductors). The cut-off energy is E c , the attractive interaction 
is V < 0, and the normal density of states per spin at the Fermi energy is denoted by 
Ao . Impurity scattering is treated in Born approximation which leads to the following 
self-energy: 

t(E,r) = ^(g(E,p F ;r)) pp (7) 

with the scattering time t. 

It is important to note that in principle the quasi-classical theory is only valid if all 
energy scales are small compared to the Fermi energy, E F ; for most superconductors 
this is the case as T c <C E F . However, for many ferromagnetic materials the exchange 
energy is of the same order of magnitude as E-p- Strictly speaking this theory can 
therefore only be applied for rather weak ferromagnets. Strong ferromagnets have 
been treated in some special cases: superconductors in proximity to half-metals 
(where only one spin-channel is metallic) or m contact to strong ferromagnetic 
insulators |201 IS] have been examined by extensions of the quasi-classical theory. 



02 {■n(p' F )g{E,p F ;r)) 



PF 



(6) 



3. Results 



First we investigate the oscillations of the pairing function in a quite simple system 
which we can treat analytically: we consider a ferromagnetic layer of thickness dp 
attached to a bulk superconductor (dg — * oo, see figure ^) without disorder; for 
simplicity we consider the case of a spherical Fermi surface in both materials, and 
an identical Fermi velocity, v F . Furthermore we assume a completely transparent 
interface at x — 0, a specular surface at x — dp, and a spatially constant order 
parameter in the superconductor, A(pp,r) = A(pp). Then the normal part of the 
Green's function in the ferromagnetic metal reads 

, F „ , _ 1~ e^n( E -^a{E,pp)P{E,p> ¥ ) 
3TT/U(£; ' PF) ~ l + e^m^^p^i^) (8) 

with 



^ E ^ = a(pV) ' (10) 




Figure 2. The density of states at E = is shown as a function of the 
ferromagnetic layer thickness with £p = irvp/I. The s-wave and the unrotated 
d-wave case are almost identical with maxima at dp = (1/2 + fc)£p/2, where 
k is an integer; for a rotated d-wave order parameter the maxima are shifted 
to dp = fc£p/2. For the plot we added a finite imaginary part to the energy 
(E — * E + i0.0lT c ). 



and 

vi/i(EiP F ) = ; 11 

\(Pf) — ^fPf/p^f {pp~- Fermi momentum parallel to the y-z-plane) is the length of 
the classical trajectory in the ferromagnetic layer. Note that p' F is uniquely determined 
by pf for a specular surface since the parallel momentum is conserved, p F = p' F ^. 

The angle-resolved density of states in the ferromagnetic layer can be expressed 
in terms of the normal part of the Green's function, 

M(E,p F ) = ^AfoEe[g n {E,p F )+ gil (E, PF )}. (12) 

We now consider the angle-averaged density of states at E = as a function of dp , 
which is normalised by the ferromagnetic length, £p = ttvf/I. In figure [21 we present 
the results for the s-wave case (A(p F ) = 1.76T c0 ) as well as for the unrotated (A(p F ) = 
2.14T c0 ?7 ( j(PF)) and for the 45°-rotated d-wave cases (A(pp) = 2.14T c0 ^(pf)); we 
choose an exchange field of / = lOT^o- 

The zero-energy density of states for an s-wave and an unrotated d-wave order 
parameter behaves quite similar up to minor deviations: we find maxima at gJf = 
(1/2 + fc)^F/2, k S No. At these values of d F Andreev bound states with a large 
spectral weight exist in the gap region; i.e. in the gap region the density of states 
is enhanced compared to the normal state value. For a 45°-rotated e?-wave order 
parameter we find maxima at g?f = The reason for this shift is that the 

quasi-particles acquire an additional phase due to the scattering at the surface which 
changes the sign of the order parameter. In particular for = the commonly known 
zero-energy bound states occur at surfaces [2*2*j . 

The density of states has been studied before El 1231 including surface 
roughness and a finite transparency of the superconductor-ferromagnet interface. The 




Figure 3. The critical temperature of a thin superconducting layer on top of a 
bulk fcrromagnet (I = 10T c oOi <^F ~~ * °°) as a function of the layer thickness, e£g, 
which is normalised by the superconducting coherence length £s (upper panel). 
For comparison the critical temperature is also shown for / = (lower panel). 



oscillatory behaviour has also been observed experimentally |141 ITS] . A discussion of 
the non-magnetic case can be found in |24|. 

In the following we focus on the critical temperature of bi-layer systems. As 
before we assume a completely transparent interface between the ferromagnet and 
the superconductor; furthermore the Fermi surfaces are supposed to be identical in 
both materials. The ferromagnet is described by the exchange energy, /, and the 
impurity scattering strength, l/2rp; we choose these parameters to be I = 10T Cj o and 
l/2rp = 5T Cj o- This is a reasonable choice having in mind ferromagnetic metals like 
Fe or Ni which are well-described by the relations T c -C /, and 1/2tf -C I- 

First we study the critical temperature of the system with an s-wave 
superconductor, which we assume to be dirty, l/2rg = 10T Cj o- The coherence length 
at zero temperature, £s, is given by £s = V^oh ~ 0.53£o where £o = vf/Ao7t (Ao = 
1.768T c .o) is the BCS-value of the coherence length for a clean superconductor, and 
— t~svf is the mean free path. In figure|!]]we present the critical temperature, T c , for 
a superconducting layer on a bulk ferromagnet (dp — > oo) as a function of its thickness, 
ds- The critical temperature decreases with the thickness of the superconductor, and 
below a critical layer thickness, ds = 1.37£s, the superconductivity vanishes. Now 
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Figure 4. The critical temperature of a bi-layer with fixed dg as a function of the 
ferromagnetic layer thickness dp (upper panel). For all order parameter types T c 
is oscillating. The parameter d$ is chosen such that T c = 0.54T Ci o for dp — * co; 
i.e. dg equals 1.72£g,1.22£sj ar >d 5.52£g for the s-wave, unrotated and rotated 
d-wave situations, respectively. For comparison the critical temperature is also 
shown for the unrotated case, 1 = 0, where no oscillations occur (lower panel). 



we fix ds = 1.72£g, for which T c = 0.54T c ,o when dp — > oo; for this thickness of 
the superconducting layer we examine the critical temperature as a function of the 
ferromagnetic layer thickness, dp. We find an oscillation of the critical temperature 
when varying the thickness of the ferromagnet (see figure 

The oscillations can be explained as follows: since no current can flow across the 
surface at x = dp, the pairing function has to obey particular boundary conditions, i.e. 
its derivative normal to the surface has to vanish. In the presence of a spin exchange 
field the pairing is spatially oscillating with wave length £p = Trvp/I [7]. At the 
maxima of T c the thickness of the ferromagnetic layer, dp, is such that the boundary 
conditions are fulfilled quite naturally, whereas minima occur if the pairing function 
has to be suppressed considerably to fulfil the boundary conditions. Therefore the 
distance of two neighbouring minima is expected to be of the order of 0.5£f, which 
is the same periodicity as observed before for the density of states (see figure 0). In 
our numerical calculation the first two minima of T c can be found at dp = 0.16£f and 
dp = 0.62£f (s-wave case, see inset in figure^J. Their distance is 0.46£f which is close 
to the expected value. These results are consistent with other theoretical findings OE] 



Superconductor- Ferromagnet Bi-Layers 
1.04 



i | i | i | i | i |- 

— d F = 0.16£ F 
---- d F = 0.39£ F 

- - d F = 0.62£ F _ 




0.98 1 ' 1 ' 1 ' 1 ' 1 ' u 

0.25 0.5 0.75 1 1.25 

E/A 



Figure 5. The density of states at x = dp for an s-wave symmetry of 
the order parameter. These values of dp correspond to the first two minima 
(dp = 0.16£p, 0.62§f) and the first maximum (dp = 0.39§f) of the critical 
temperature. 




Figure 6. The density of states at x = dp for an unrotated d-wave symmetry of 
the order parameter. These values of dp correspond to the first two minima 
(dp = 0.20,0.67§p) and the first maximum (dp = 0.43£p) of the critical 
temperature. 



which could also be fitted to experimental observations ■ For comparison we also 
present the critical temperature for a non-magnetic metal layer (1 = 0); as expected 
superconductivity is suppressed less effectively without magnetism (see figure yj. For 
this case, no oscillations with the metal layer thickness are observed (see figure 0}, 
which is an obvious result considering the discussion in relation with (|^|)- (|ll|l . 

As previously discussed, the oscillating behaviour of the pairing function can also 
be observed in the local density of states at x = d F . We calculate the density of states 
for those values of d F where the critical temperature has a maximum or a minimum. 
The results are presented in figure [5J For the minima of T c (at d F = 0.16,0.62£p) 




Figure 7. The density of states at x = dp for a 45°-rotated d-wave symmetry 
of the order parameter. The left axis applies to the density of states at the 
first maximum of T c (dp = 0.15£p), whereas the right axis applies for the first 
minimum (dp = 0.39£f) and the second maximum (dp = 0.62£p). 



the density of states in the gap region is enhanced compared to the normal state; 
this is related to sub-gap Andreev bound states which exist inside the ferromagnet. 
At the maximum of T c (dp = 0.39£f) we find a gap-like structure of the density 
of states; i.e. in the gap region the density of states is smaller than the normal 
state value. Altogether we find that bound states in the gap region tend to suppress 
superconductivity, and can be related to minima of T c . 

Now we will turn to the case of d-wave superconductors. For a d-wave symmetry 
of the order parameter the bulk value of the critical temperature, T CjU , is suppressed 
by non-magnetic impurities, and T c is given by |25| 



4nT Ct0 T S t 

where T c .oo is the critical temperature for a clean sample and ^(x) is the digamma 
function. In the following we choose the impurity scattering inside the superconductor 
to be l/2rs = 0.1T Cj oo which leads to T Cj o = 0.92T Cj oo; the bulk order parameter at 
zero temperature is Ao = 2.02T Cj oo which is smaller than its value in the clean case 
(A o = 2.14T C) oo). As the disorder is small the superconducting coherence length is 
given by £ s » £ = u F /A 7r. 

It is well-known that the behaviour of d-wave superconductors at boundaries 
depends crucially on the orientation of the order parameter with respect to the 
boundary. Therefore we compare the case where the order parameter is rotated by 
45° with the unrotated case. 

An unrotated d-wave superconductor is expected to behave similar to an s- 
wave superconductor. The reason is that along the classical trajectories the order 
parameter does not change its phase. And indeed the behaviour of an unrotated d-wave 
superconductor and an s-wave superconductor is quite similar: the superconductivity 
of a thin layer on a bulk ferromagnet is suppressed completely when its thickness is 
below ds — 1.22£s (see figure O ■ If we fix the superconducting layer thickness to 
ds = 1.94£s (which leads to T c — 0.54T Cj q for dp — * oo) and vary dp we find an 
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oscillating behaviour of T c as before. The first two minima are at dp — 0.20£p and 
dp = 0.67£p, and the first maximum can be found at dp = 0.43£p (see figure @J; 
the distance between the first two minima is 0.47£p- Of course some quantitative 
differences exist which are mainly due to the nodes of the d-wave order parameter. 
The critical temperature for a non-magnetic metal layer (I = 0) behaves similar to 
the s-wave case: for an infinite metal layer T c lies above the value for the magnetic 
case (see figure 0J; the oscillations with a varying thickness of the non-magnetic metal 
vanish (see figure 0J. 

For those values of dp which are related to a maximum or a minimum of T c the 
density of states shows qualitatively the same behaviour as for the s-wave case; i.e. 
for small energies the density of states at minima of T c is enhanced compared to the 
normal state value whereas a gap-like structure can be observed at the maxima of T c 
(see figure EJ). 

The situation changes drastically if the d-wave order parameter is rotated by 
45°: Surfaces are now pair-breaking as the quasi-particles are scattered so that a 
sign change of the order parameter occurs along their trajectories |22| : and it is well- 
known that at specular surfaces the order parameter vanishes. As a consequence 
the suppression of superconductivity is much stronger than for the cases discussed 
before. For a superconducting layer which is on top of a bulk ferromagnet this 
can be seen in figure [3J in particular the critical thickness of the superconductor, 
below which superconductivity vanishes, is ds — 5.18£s, which is much larger than 
in the previous situations. We now fix the layer thickness of the superconductor to 
ds = 5. 52£s so that T c = 0.54T c ,o for dp — ► oo. When analysing T c as a function 
of dp we find a completely different behaviour than before: for dp = the order 
parameter at the pair-breaking (specular) surfaces, x = —ds and x = 0, must be 
zero. This suppression of superconductivity leads to a vanishing critical temperature 
of the bi-layer for dp < 0.06£f- The pair-breaking at x = is weakened when the 
ferromagnetic layer thickness increases, and for dp > 0.06£f the critical temperature 
becomes finite. When further increasing dp the critical temperature is oscillating as 
in the previous cases, but now starting with a maximum at dp = 0.15£f; the first 
minimum can be found at dp — 0.39£f, and the second maximum at dp = 0.62£p- 
The difference between the first two maxima is 0.47£f as for the unrotated order 
parameter. It is not surprising that the periodicity is not affected by the rotation 
of the order parameter because the oscillation of the pairing function is an exclusive 
result of the ferromagnetic exchange energy 7, and is independent of details inside 
the superconductor. For an infinite non-magnetic metal layer (I = 0) the critical 
temperature of the rotated d-wave superconductor remains unchanged (see figure [21), 
as the suppression of superconductivity is not dominated by the metal layer but by the 
surface of the superconductor at x = —ds as discussed above. The critical temperature 
shows no oscillations with the thickness of the non-magnetic layer (see figure 0} . 

The density of states at x — dp shows a clearly different behaviour. In particular 
for dp = zero-energy Andreev bound states exist at the surfaces due to the sign 
change of the order parameter for scattered quasi-particles. For the first maximum 
of T c (dp — 0.15£f) the density of states is suppressed below the normal state value 
(see figured) but a remainder of the zero-energy bound state is still observable. For 
the first minimum of T c (dp = 0.39£f) the density of states for small energies (E w 0) 
is enhanced compared to the normal state value, and for the second maximum of T c 
(dp = 0.62£f) it is suppressed, which, however, can hardly be seen in figured It is 
remarkable that the density of states in the gap region has more structure here than in 
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the previous case. This is due to the strong angular dependence of the order parameter 
close to those directions which are perpendicular to the interface. Altogether, for an 
45°-rotated d-w&ve order parameter, we also find that the minima of T c are related to 
an enhanced density of states in the sub-gap region, and vice versa for the maxima of 
T 

4. Conclusion 

We have studied superconductor-ferromagnet bi-layers for s-wave and <i-wave 
superconductors. In all discussed cases we observed an oscillating behaviour of the 
density of states as well as of the critical temperature when varying the thickness of 
the ferromagnetic layer. The origin of these oscillations is the exchange field in the 
ferromagnetic metal which leads to a Larkin-Ovchinnikov-Fulde-Ferrell-like state with 
a spatially oscillating pairing function. In particular we find that the density of states 
in the ferromagnetic layer is enhanced in the gap region when its thickness leads to a 
minimum of the critical temperature. When the critical temperature has a maximum 
the density of states in the ferromagnetic layer has a gap-like structure. 

Comparing the different order parameter symmetries, we observe a similar 
behaviour for the s-wave and the unrotated d-w&ve cases. The critical temperature 
as a function of the ferromagnetic layer thickness, dp, decreases for small values of dp 
(dp < 0.2£p) and shows oscillations around an asymptotic value of T c when further 
increasing dp. In the s-wave case these findings are in agreement with previous 
theoretical studies [HI Ej , and are also experimentally confirmed 0] |5] . 

This behaviour is considerably modified if the ci-wave superconductor is rotated 
by 45° with respect to the surface. In this case superconductivity may even vanish 
for very thin ferromagnetic layers. If c?f exceeds a critical value, superconductivity 
can be restored and oscillations around the asymptotic value of T c are observed. 
This difference in behaviour is due to the sign change of the order parameter for 
quasi-particles which are scattered at the surface, which leads to a suppression of 
superconductivity. It would be most interesting to check our results for e?-wave 
superconductors also experimentally. 
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